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Abstract 

The dually transformed path integral of four-dimensional U{1) lattice gauge theory is 
used for the calculation of expectation values in the presence of external charges. Applying 
the dual simulation to flux tubes for charge distances up to around 20 lattice spacings, we 
find a deviation from the behaviour of a dual type-II superconductor in the London limit at 
large distances. The roughening of the flux tube agrees with the effective string picture of 
confinement. Further we show that finite temperature effects are negligible for time extents 
larger than the charge distance. Finally we analyze the different contributions to the total 
energy of the electromagnetic field. The parallel component of the magnetic field turns out 
to have the strongest influence on the deviation from the linear behaviour at small charge 
distances. 

1 Introduction 

The mechanism of quark confinement within hadrons is still an open problem. Numerical simu- 
lations of QCD on space-time lattices demonstrate the confinement phenomenon but could not 
yet clarify completely the mechanism for the formation of a gluonic flux tube. Effective theo- 
ries have been created to improve our understanding of this mechanism. An intuitive physical 
picture is the dual superconductor model Dynamically generated colour magnetic monopole 
currents form a solenoid which squeezes the chromoelectric flux between quarks into a narrow 
flux tube. The energy of this flux tube increases linearly with its length and confinement is 
achieved. Another approach is the effective string theory |§, ^ which should describe flux tubes 
at large quark distances. It is an interesting goal of lattice gauge theories to test the validity of 
such effective models and therefore deepen our understanding of the confinement mechanism. 

U{1) lattice gauge theory in four dimensions undergoes a phase transition at strong coupling 
and has been widely used as a prototype of a confining theory. In U{1) theory the condensation 
of magnetic monopoles is made responsible for the confinement of electric charges. The results 
of lattice simulations have been interpreted in terms of the dual superconductor picture of 
confinement, e.g. the validity of a dual London relation has been checked Q, and the U{1) 
results for electric fields and monopole currents were found to agree with the predictions of a 
classical effective model considering also string fluctuations |^. Most of the obtained results 
are expected to hold also in non-Abelian gauge theories, at least qualitatively The data 
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accumulated so far, however, for larger flux tube lengths are not of sufficient accuracy to perform 
a quantitative analysis in terms of effective models for the formation of a gluonic flux tube. 

Furthermore, it was realized many years ago ||^] that one can perform a duality transformation 
of the path integral in compact Abelian gauge theories. In this way a new partition function 
is obtained which can be regarded as a limit of the dual non-compact Abelian Higgs model Q. 
We will use the duality relation for expectation values in the presence of external charges to 
clarify the connection between dually transformed U{1) theory and the dual superconductor 
model. Besides this support for the interpretation of the results, the dual theory provides a very 
efficient tool for simulating flux tubes. This fact has already been realized in refs. but to 

our knowledge no dual simulations of four-dimensional U{1) lattice gauge theory with external 
sources have been performed yet. 

The organisation of this article is as follows: In section 2 we review the duality transformation 
of U{1) lattice gauge theory for the Wilson action with special consideration of fields and their 
squares in the presence of external charges. The consequences for the interpretation in the dual 
superconductor picture are discussed by connecting the results to a dual Higgs model in section 
3. After pointing out the advantages of dual simulations of flux tubes in section 4, we present our 
results for electric field strength and magnetic current in the presence of a charge-anticharge- 
pair in section 5, analyzing the flux tube as a function of spatial lattice size and temperature. 
In section 6 we return to the dual superconductor picture, subjecting our data to a quantitative 
comparison with the London model. In section 7 the dependence of the flux tube width on its 
length is investigated and compared with the prediction of an effective string theory. Finally 
in section 8 we calculate the total electromagnetic energy and analyze the various components 
contributing to it. 



2 The duality transformation of U{1) lattice gauge theory 

We begin with the usual prescription for simulating U{1) lattice gauge theory. Expectation val- 
ues of physical observables O in the presence of a static charge pair at distance d are determined 
by the correlation function 



( L(O)L+(d)0(x )) 
(L(0)L+(ci)) 



where L(r) is the Polyakov loop L{r) = f^x=(r,fca),^=o ^^'^ the angle brackets denote the 

evaluation of the path integral using the standard Wilson action 

Sw = P [I - cos{deU^,]. (2) 

X,fl<U 

(3=l/e^ is the inverse coupling, and {d9)x^^u is the discretized exterior derivative of the phases 
9 of the link variables Ux,f_i = e*^""-'' and is assigned to a plaquette of size a^. Below we will 
use the notation of lattice differential forms and therefore suppress indices. The observables O 
we are interested in are the electric field strength, the magnetic current and the electric and 
magnetic energy densities. For the field strength F we use the identification a^eF = sm{d9) 
which was shown to fulfil the electric Gauss law for the Wilson action . The magnetic current 
can be constructed via the dual Maxwell equations. The squared fields are calculated in the 
usual way by means of the operator (5cos{d9). 
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, (4) 

Sk=0 



Let us briefly review the most important steps of the duahty transformation. We start from 
the U{1) partition function using the Wilson action which reads [|l3| 

Z = 11 r de e^p{-f3 cos dO)}, (3) 

J —IT 

where the product is over all links on the lattice, and the sum in the exponent is over the 
plaquettes. Expanding the Boltzmann factor for each plaquette in a Fourier series one can write 

Z = l[r de HE exp{i(A;,d^)}e-%||(/3) = (27r)4^nE ^'^hkwiP) 
k k 

where we have used the plaquette variables k and the property of the inner product (fc, de) = 
{5k, e), and then performed the integration over the links e. N is the number of lattice sites, and 
the product has to be performed over all plaquettes. k is an integer valued 2-form corresponding 
to the field strength F. The constraint 5k = describes the absence of electric charges, 5F = 
Jg = 0. We solve this constraint by introducing an integer valued 3-form I with 

k = 5l, (5) 

z=(2^)4^nE ^-'hsiwiP) = (2vrrnE ^-'imm^ (6) 

where we have finally switched to the dual link variables representing an integer valued dual 
potential. The product in the partition function extends over the links of the dual lattice. 

The next step is the inclusion of external charges into the U{\) path integral. and £_ 
denote the world-lines of a static charge pair. Starting from 

ZQQ = ]ir de eM-f3T.(^-cosde)}l[e-^'l[e^', (7) 

one arrives at an expression equivalent to with a modified constraint 5k = ^1 for links 
belonging to C± (equivalent to e5k = Jg). This constraint is solved by an additional integer 
2-form n: 

k = 51 + n, with e5n = Je, resp. *k = d*l + *n, with ed*n = *Je. (8) 

External electric currents are represented by the boundaries of a dual Dirac sheet in the dual 
theory. The dually transformed path integral therefore reads 

Zqq = (27r)4^ nE e-^h<iH+*nm, (9) 

where the product is performed over all links of the dual lattice. The dual Dirac sheet *n can 
be fixed as the minimal straight area between the two Polyakov loops, while the integer gauge 
field *l is kept as dynamical variable and updated in simulations. Fluctuations in *l result in 
fluctuations of the flux string d*l + *n. 

Further, we will perform the duality transformation for the correlation function (||). For 
the determination of the flux Foq through a given plaquette Do, we have to insert the operator 
V^sind^Dp into the path integral in addition to the Polyakov loops. So we can write 

< >QQ = z-^h r d0 HE n^-"n^" # fn-" -nH >< 

X exp{i{5k,e)} hk\\{/3) . (10) 
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Taking the imaginary part reflects the change of sign of the field strength under an exchange 
of charge and anticharge, i. e. under complex conjugation. Our constraint for 5k is further 
modified: 6k has to be non-zero along the electric current and the boundary of the considered 
plaquette Oq. Solving this constraint by a redefinition of k on this plaquette (/c — > /c it 1) and 
using the identity for modified Bessel functions 

kim = ^[ik^m-h+im, (n) 

one obtains 

<Fa,>QQ = Z^l (2vr)4^ [IE ;^*^°o ^''^ hdH+*n\\m, (12) 

where *k = d*l + *n as before. Note that in this dual expression the validity of Gauss' law for 
electric charges is checked easily because of d*k/y/P = *Je- 

Performing the duality transformation for the expectation value of /? cos dOn^ in an analogous 
way leads to 

<<>Q0 = {2^r IIE ^ff^f e-^Im+^m - <F'>.ac. (13) 

We would like to mention that an identity equivalent to ( p^ ) was already derived in ref. Q for 
calculating the electromagnetic energy density in three-dimensional U{1) lattice gauge theory, 
using the polymer formulation of the path integral. 



3 A comparison to a dual Higgs model 

The results of C/(l) simulations in the confinement phase verify the dual superconductor picture 
0, ^: The electric flux between a charge pair is squeezed into a flux tube, encircled by monopole 
currents acting like a coil. We want to discuss how this scenario can be described by the dual 
non-compact Abelian Higgs model, in order to establish an exact connection to the dually 
transformed lattice U{1) theory presented above. The dual Higgs model is written in terms 
of dual potentials *A assigned to links on the dual lattice (the *-notation is again used for all 
fields which are located on the dual lattice) and of a dual Higgs field which is located on 
the dual sites. Electric charges in this dual model correspond to monopoles (assigned to three- 
dimensional cubes on the dual lattice) and can be introduced by the dual of a Dirac string in 
the same way as in the last section. On the original lattice the world lines of electric charges 
define the boundary of a continous surface (the Dirac sheet): 

{1 for plaquettes on the surface , . 

elsewhere 

The field strength *F is given by g*F = g d* A + 27r*n, where d*A is the exterior derivative of 
the 1-form * A, and g = lixje is the magnetic coupling. In the London limit the Higgs field *$ 
is constrained to *$ = exp(z*x) and therefore the coherence length vanishes. Including electric 
sources the corresponding action reads 

S = /3 G(5dM + 27r*n) - 7 ^ cos((i*x-5*^), (15) 

plaquettes links 



4 



where /? = l/g^- The function G{g*F) determines the non-compact gauge field action and 
approaches G{g*F) — > g'^ *F'^/2 in the continuum hmit. The second part of the action is the 
compact lattice version of the gauge invariant kinetic term of the Higgs field The ratio /3/7 
has the meaning of the bare squared London penetration length A^. The Higgs current *Jh, 
given by 

g*JH = l {[d*x]-g*A), (16) 

is constructed in such a way that the fluxoid quantisation *F + X^d*JH = = *me Q is 

valid for every single field configuration. The integer two-form *m characterizes the physical 
fluxoid string which is allowed to fluctuate in simulations due to the compactness of the phase 



X of the Higgs field. [d*x] in eq. (16) indicates the reduction of d*x to the interval (— 7r,7r]. 
Therefore (i[(i*x] may be non-zero for topologically non-trivial field configurations, and 27r*m = 
27r*n + d[d*x]- If the value of [d*x] on one dual link flips, i. e. the phase difference between 
neighbouring points exceeds the range (— 7r,7r], this effects six dual plaquettes corresponding to 
a closed cube on the original lattice. If such a cube has one plaquette in common with the Dirac 
sheet (n = 1), 27r*n and d[d*x] may cancel each other, changing the shape of the physical fluxoid 
string. 

This Higgs model can be regarded as four-dimensional generalisation of the static classical 
effective model investigated in ref. B. It can also serve to extend numerically the predictions 



from the dual QCD ansatz |11] or from similar analytical calculations in effective models [12|, 
especially in investigating fluctuations of the fluxoid string. For the U{1) gauge theory, it is 
even possible to establish an exact connection at the level of the path integral. Let us assume 
7 — > oo which imposes a new constraint on the fields: Choosing a gauge where the phase x of 
the Higgs field in eq. (|l5|) is zero, we realize that the integral over *A for each dual link reduces 
to a sum over an integer */, according to the constraint *A = 27r*/. Up to a /3-dependent factor 
the 7 — > oo limit of the partition function of this Higgs model equals the dually transformed 
path integral (^ of U{1) gauge theory, including two Polyakov loops. This duality property is 
well-known for the Villain form of the action [Q. In the case of the Wilson action the function 
G{g *F) describing the gauge field action in the corresponding Higgs model is determined by the 
modified Boltzmann factor in eq. (Q) 



exp 



-PG{g*F) =e-^/||„,^||(/?), (17) 



and ^= l/(47r2/3). 

Especially interesting for us is the interpretation of expectation values in ^7(1) in terms of 
the dual Higgs model. Let us first regard the field strength in U{1) derived in eq. (|12|): It also 
agrees with the field stength of the dual Higgs model in the 7 — > cxd limit which reads 

g*F 2Trd*l + 2-K*n ^ *F ^ e{d*l + *n) = -^*k. (18) 

This means that the field strength expectation values and due to the validity of the extended 
Maxwell equations also the magnetic currents in U{1) gauge theory agree with the corresponding 
fields and currents in the dual Higgs model in the limit of both zero coherence length and zero 
penetration length. From this point of view a ^7(1) fiux tube looks rather like an infinitely thin 
fiuctuating string than like a fiux tube in a classical dual superconductor. 
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This correspondence to the Higgs model does not hold for energy densities. The dually 
transformed expectation value of the squared field strength, eq. (|l3[), obviously disagrees with 
the corresponding result in the 7 ^ oo limit of the dual Higgs model. This is not astonishing, 
since the expectation values of the squared fields are connected to the vacuum fluctuations 
which show the opposite behaviour in the dual theory, where the confinement phase is the 
weakly coupled phase and fluctuations become more important with increasing (3. Hence, a 
direct interpretation of the U{1) energy density in the dual Higgs model is not possible, which 
makes it difficult to distinguish by analytical considerations between dual type-I and type-II 
sup er conductivity. 



4 Simulating the dually transformed U{1) theory 

As has been pointed out above, one has to be careful when comparing U (1) lattice gauge 
theory to a dual Higgs model on a "microscopic" level. Nevertheless, simulating the dual U{1) 
gauge theory gives us the possibility to calculate any expectation value in the presence of static 
charges with significantly higher precision. This has several reasons: The confinement phase is 
the weakly coupled one in the dual theory, therefore we have less quantum fluctuations in the 
region of interest. Simulating the dually transformed path integral (|l2|) resp. ( p^ ) is also much 
faster because of the integer valued gauge fields. Most important, however, is the fact that it 
is not necessary to project the charge-anticharge state out of the vacuum: Charge pairs with 
arbitrary distance can be simulated with equal accuracy, as well as multiply charged systems 



MM- 



The updating procedure itself is very simple: As already discussed in the previous section, the 
dual Dirac string sheet *n connecting the two Polyakov lines is held fixed, while the integer gauge 
field located on the dual links, is changed by ±1 during the updates. We have to mention 
that this is essentially the same procedure as the "change-a-cube" algorithm used in ref. |p. The 
only difference for the three-dimensional model used there is that the dual potential is located 
on the dual sites instead of the dual links. In ref. 1^] the electromagnetic energy distribution 
was calculated to extract the string tension in a range of weaker couplings as total field energy 
per length. However, this quantity must not be confused with the string tension calculated from 
usual correlation functions. 

Recently, a detailed analysis of simulating dual U{1) lattice gauge theory in three dimensions 



was performed in ref. [16|. The authors found that at weak coupling the advantages are out- 
weighed by large autocorrelation times as the dual model becomes disordered. We also observed 
the increase of the autocorrelation time and correspondingly of the error for increasing values 
of /3 - with a peak at the phase transition of the four-dimensional theory. In the deconfinement 
phase there are large fluctuations, but the autocorrelation time is smaller than around the phase 
transition. Since our goal is the investigation of the conflnement mechanism, we take the main 
advantage from simulations in the confinement phase, anyhow. In this context we want to men- 
tion that we used a duality transformation exactly valid for the whole (5 range, not only in the 
weak coupling limit as in ref. [^]. Further, as stated above, the power of dual simulations lies 
in the inclusion of the sources in the action, which already has been realized in ref. IC]. 

For each charge distance d a separate simulation has to be performed, and of course for 
the subtraction of the vacuum expectation value in (|l]) a distinct run is needed. Therefore 
one gets a good estimate of the reliability of the errorbars. For smaller distances on a 8^ x 4 
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lattice an accurate comparison was performed between expectation values from ordinary and 
from dual simulations 1 14 1 . The only difference in the observed results is the effect of the periodic 
boundary conditions which is absent in the dual model: There can be no electric flux over the 
period per construction. For d > 3a the dual simulation is clearly more efficient, because in 
ordinary simulations the signal decreases exponentially with the charge distance d, while in the 
dual simulation it stays constant (the necessary amount of computer time still depends of course 
on the lattice size used). 

For the updating procedure we use a standard Metropolis algorithm, which has proven 
sufficient for the investigated range of couplings, although a nonlocal cluster algorithm would 
certainly lead to a further increase in performance. Our typical measurements are taken from 
5 • 10^ — 2 • 10^ configurations, where each 10th configuration is evaluated. For determining 
the errorbars we took blocks of 100 evaluated configurations. For the results shown in the next 
sections we used the coupling /3 = 0.96, where the system is clearly in the confinement phase. 
On the other hand the fluctuations in the dual theory are already large enough to produce a 
flux tube of the desired transverse extent. If one likes to compare with the QCD flux tube, we 
can require that our string tension equals the physical value, which gives a lattice spacing a of 
roughly 0.25/m. 



5 The profile of fiux tubes between static charges 

We present now our numerical results for the electric field and magnetic current distribution 
in the presence of a static charge pair. The great advantage compared to the results reported 
in ref. ^] is that we are now able to investigate much larger charge distances and much lower 
temperatures. 

Before discussing the temperature dependence of the fiux tube profile, let us comment on the 
necessary spatial size of the lattice. There is the simple condition that the spatial extent of the 
lattice has to be greater than the extent of the considered fiux tube both in longitudinal and in 
transverse direction in order to prevent strong finite size effects. For our simulations we used a 
spatial extent of 8^ for charge distances d = la — 5a, 12^ for simulating d = 5a — 8a, and 16^ for 
d = 8a — 12a. Further we performed some simulations on larger lattices; the largest investigated 
charge distance was d = 22a (on a 26^ x 32 lattice). Results for this distance are displayed in a 
3d-plot in fig. m. We focus on the longitudinal component of the electric field strength and on 
the azimuthal component of the magnetic current which encircles the fiux tube like a coil. It can 
already be seen in these plots that a flux tube of constant thickness is not observed even for this 
large charge distance. The increase of the flux tube width with the distance between charges 
may be explained by roughening effects and will be discussed in section |^. The right flgure 
demonstrates that the signal for the magnetic current drops in the middle between charges, its 
radial dependence in the symmetry plane exhibits a maximum at finite R. This behaviour is 
interesting for the dual superconductor interpretation as will be discussed in the next section. 

In all our simulations we represent static charges by Polyakov loops and therefore have to 
deal with finite temperature effects due to the finite time extent of the lattice. It is therefore 
important to separate the occuring thermal fluctuations from quantum fluctuations which are 
essential for the extended fiux tube. This is achieved by analyzing the dependence of the results 
on the time extent NT, as is shown in fig. ^ for the electric field profile and in fig. ^ for the 
magnetic current profile in the symmetry plane between charges for various charge distances 
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Figure 1: The flux tube between static charges for a distance of 22 lattice spacings. The left plot 
shows the parallel component of the electric field, the right plot shows the azimuthal component 
of the magnetic current, both in a half plane through the QQ-axis. Note that in the symmetry 
plane between charges the magnetic current takes its maximum value at a transverse distance 
of 1.5a, not at R = 0.5a. 

d. From the results we conclude that if the time extent NT of the lattice exceeds the charge 
distance d finite temperature effects can safely be neglected. For smaller NT, however, we realize 
a significant widening of the fiux tube, connected with a decreasing magnetic current. Prom the 
point of view of the dual simulation, this widening can easily be explained. Quantum fluctuations 
of the world sheet of the string correspond to cube-like excitations on the original lattice. At 
finite temperature there also appear time-independent excitations which wrap around the lattice 
in time direction and therefore increase the effective fiux tube width. These excitations vanish 
if NT is large compared to the length of the fiux tube. 

After the discussion of the techniques used for the simulation of long flux tubes and the 
dependence of the results on the temperature, we will now turn again to the dual superconductor 
interpretation. 

6 A closer look at the dual London equation 

In the presence of a pair of static charges, the classical solution of the Higgs model (|l^) yields 
an electric field and a magnetic current distribution according to 

rot^(r) = -Jm{r), (19) 
E{r) = X'^votJm{r) + £{r), (20) 

where the (straight) string £{'r) connects the point charges, carrying a quantum of electric fiux 

e. Integrating eq. (|20|), a generalisation of the London equation, over a surface intersecting the 
string gives the fiuxoid quantisation. For an infinitely long string eq. (|l^) and eq. ( pO| ) yield an 
electric field profile of the form E{R) = {e/2TT\'^)Ko{R/ X). 
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E(d=U) E(d=5a) 




Figure 2: The temperature dependence of the electric field profile in the symmetry plane between 
charges for various charge distances d: If the time extent NT of the lattice (corresponding to the 
inverse temperature) is smaller than the charge distance, the width of the flux tube is strongly 
temperature dependent. 

Jm(d = 3a) Jmid = 12a) 




Figure 3: The temperature dependence of the magnetic current profile in the symmetry plane 
between charges for charge distances d = 3a and d = 12a: If NT is smaller than d, the signal 
for the current gets significantly lower which means a widening of the flux tube. For d = 12a 
the maximum current is observed at the transverse distance R = 1.5a, like in the d = 22a case 
shown in fig. |T]. 
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E\\(d = U) E\\(d=12a) 




Figure 4: Comparison of the U{1) electric field profile with the numerical solution of Maxwell 
and London equation for a straight fluxoid string. The free parameter A is fitted. For small 
charge distances (d = 3a, left figure, Xfu = 0.75a) good agreement can be achieved; for larger 
distances {d = 12a, right figure, Xfn = 1.50a) it is not possible to describe the U{1) flux tube 
behaviour. 

These classical equations were used in ref. for comparison with the results of [/(I) lattice 
simulations. The free parameter A (the London penetration length) could be fitted by comparing 
E and rot off axis, and also the fluxoid quantisation was verified. An extension of this analysis 
was performed in ref. by including fluctuations of the string £{r) in the above model. Further 
the effective model was solved numerically on a lattice of same size, to avoid being misled by 
lattice artifacts. Due to the exponential falloff of the signal, however, the U{1) data analyzed 
so far contain only charge distances up to d = 3a. We will now consider the data obtained in 
the dual simulations, in order to see if the above interpretation holds for longer flux tubes. 

In the classical model with a straight string £{r) one would expect that 

• the flux tube width as a function of d soon approaches a constant, this means that the 
fitted parameter A should not depend on d, 

• A^ is constant everywhere off axis, 

• the azimuthal magnetic current encircling the axis should exhibit a maximum in the sym- 
metry plane between the charges, and not near the charges. The radial dependence should 
exhibit an exponential falloff. 

It was shown that these requirements are not even fulfilled for a charge distance d = 3a Q, 
nevertheless the flux tube profile in the symmetry plane may be described well by the effective 
model with a static string Q. Considering larger charge distances confirms the discrepancy to 
the items mentioned above. In a comparison of electric field profiles for d = 3a and d = 12a 
(see fig. ^) it can be seen that the fitted values for A not only differ very much, but it is already 
impossible to describe the field distribution for larger charge distances by the numerical solution 
of coupled Maxwell and London equation. 

For the magnetic current distribution we already see in fig. |l| and in fig. ^ that its behaviour 
cannot be explained in the above effective model. Fig. ^ shows the corresponding classical 
solutions, again for d = 3a and d = 12a. For larger distances, there is not even qualitative 
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Jm(d = 3a) Jmid = 12a) 




1 2 3 1 2 3 

R[a] R[a] 

Figure 5: Comparison of the U{1) magnetic current profile with the numerical solution of 
Maxwell and London equation for a straight fluxoid string (A- values chosen as in fig. ^). The 
agreement is reasonable for small charge distances (d = 3o, left figure), for large distances the 
behaviour of the current distribution cannot be described by the solution of the London model 
{d = 12a, right figure). 



agreement: The azimuthal component of the magnetic current as a function of the transverse 
distance first increases and takes its maximum value at R = 1.5a. 

This phenomenon together with the other described observations may be explained within the 
dual superconductor picture by two different approaches: First, it can be regarded as evidence 
for the importance of string fiuctuations. This way was chosen in ref. (5|, where the intrinsic 
thickness of the flux tube A turned out to be much smaller than the effective flux tube width 
which increases with charge distance. The extreme scenario (A — > 0) leads us again to the exact 
correspondence of U{1) to the 7 — > oo limit of a dual Higgs model. 

The second possibility for arguing is to stay in the classical model of a dual superconductor, 
but leave the London (extreme type II) limit described so far. A strong suppression of the Higgs 
condensate in the region of electric fields may also lead to the behaviour of the current observed 
in fig. ^. A similar current distribution measured in Abelian projection of SU{2) was interpreted 
in this way in ref. |17]. Hence, one may conclude that U{1) lattice gauge theory in this respect 
behaves like a dual type-I superconductor. Stronger evidence for such a statement is expected of 
course from an investigation of the interaction between flux tubes (which also can be performed 
very efficiently in a dual simulation). Indeed preliminary results [|l^ show an attraction between 
U{1) flux tubes for the regarded value of (3. 

We may conclude that interpreting the flux tube results in terms of the dual superconductor 
picture still makes sense, but U (1) lattice gauge theory certainly cannot be regarded as a classical 
dual superconductor in the London limit if one examines larger charge distances. On the other 
hand the effective string description of flux tubes represents an independent feature of the 
lattice gauge system, and from the point of view of dual simulations this looks somehow easier 
to argue. We will use this approach in the next section for a closer look at the fiux tube width 
in dependence of the charge distance. 
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7 The flux tube width and the effective string description 



In ref. |^ the string model approximation to Wilson loop correlation functions was used to study 
the behaviour of the flux tube width for large charge distances. This model assumes a very thin 
bare flux tube subjected to quantum mechanical fluctuations. A convenient measure for the 
transverse extent of the flux tube is the squared width 

where e represents a quantity characteristic of the flux tube, for example the (chromo-) electric 
field energy density. It was found in ref. [^j that the squared width w"^ diverges logarithmically 
for d — > oo, 

w"^ = wl \n{d/dc) (22) 

with some constants Wq and dc- A numerical test of this hypothesis could not be performed up to 
now for U{1) lattice gauge theory. In the three-dimensional Z2 gauge model this prediction was 
verified in ref. ||lO| , where the authors also exploited the duality relation to a three-dimensional 
Ising model and combined their data with high precision data on the interface physics of the 
three-dimensional Ising model. They also achieved agreement with the prediction 

(where a is the string tension) for the two-dimensional free Gaussian model. There is the 
conjecture that this behaviour of the flux tube width should be observed in all gauge theories. 

It is therefore interesting to test this prediction in dual U{\) simulations where large flux 
tubes become accessible. We determine the flux tube width by means of the parallel component 
of the electric field strength E\\. Our defintion on the lattice reads: 

2 Y.iRlE\\{Ri) 

This slightly reduces the estimated error compared to the definition via the electric field energy, 
in all other aspects the results roughly agree. The sum extends over a transverse plane of 
the lattice with Ri < Rmax- For d < 12a we evaluate the sum until a transverse distance 
Rmax = 4a, for larger distances we use Rmax = 5a. In fig. ^ one can see again the influence 
of finite temperature effects on the flux tube width. The calculation with time extent 16 can 
be taken as acceptable approximation to the zero temperature case for the regarded charge 
distances and is analyzed according to the prediction ( |2^ ) in fig. ^ in analogy to ref. ||l^ . The 
logarithmic increase of the flux tube width is realized very well for U{1) lattice gauge theory. 
The constant Wq, however, disagrees with expr. ( p3|) by roughly a factor 2. 

The above data demonstrate that the thickness of U (1) flux tubes does not saturate for 
large charge distances, as expected in a naive formulation of the dual superconductor picture. 
It seems to diverge logarithmically, as predicted by the effective string description. As can be 
seen from fig. this behaviour sets in already at d = 2a, which corresponds to a physical length 
of around 0.5/m in agreement with ref. flC 
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Figure 6: The thickness of the flux tube we (calculated from the profile of the electric field) 
in the symmetry plane between charges as a function of the charge distance d for various time 
extents NT of the lattice. If d exceeds the time extent NT, finite temperature effects become 
important and the flux tube significantly widens. 
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8 Energy from the dual simulation 



In section 2 we have discussed the duahty transformation for the expectation value of /3 cos dOn^ 
in the presence of two Polyakov loops - see eq. (13) - which gives the squared field strength on a 
specified plaquette Oq. This allows for the determination of the Euclidean action density E'^+B^ 
and energy density E"^ — B^. As has been pointed out, the squared field strength calculated 
in the dual simulation differs from the squared field strength of the 7 ^ oo Higgs model. The 
contributions of all plaquettes can be summed up and give the total action respective energy 
of the system. In ordinary lattice simulations this is a very difficult task, since the difference 
between the squared field strength of the pure U{\) vacuum and the charge-anticharge state 
is small compared to the vacuum expectation value itself. Therefore it cannot be resolved for 
larger charge distances. Further, concerning the energy density there are large cancellations 
between electric and magnetic field contributions. 

In the dual simulation of U{\) lattice gauge theory it is possible to determine the total 
energy of the electromagnetic fields as a function of the charge distance. The influence of finite 
temperature effects is analyzed in fig. ^. The thermal fiuctuations increase the total energy. 
This is again due to the periodic excitations in time direction which contribute to the electric 
part of the energy only. 

The total field energy shown in fig. ^ qualitatively behaves like the free energy for a charge 
pair (as determined in usual lattice simulations), i.e. up to a constant like a Coulomb and a linear 
term. However, as already realized in ref. a quantitative comparison shows that the integral 
over the energy density differs from the potential. We mention here that it is also possible to 
determine the free energy via dual simulations; this opens the opportunity for the application 
of lattice sum rules |jl^. Such calculations are in progress and will be reported elsewhere. 

It is also interesting to look at the individual contributions to the total field energy. For 
/3 — > there are no string fluctuations, and therefore only the parallel component of the electric 
field contributes. If we consider single cube-like excitations in transverse directions, we get also 
contributions from transverse electric and transverse magnetic fields. In this case they cancel 
each other exactly and therefore do not change the total energy. For stronger fiuctuations, 
i. e. increasing roughening of the string, these components differ, and also the parallel magnetic 
component gets nonzero. The results for /3 = 0.96 are plotted in fig. |^. One can argue that 
for charge distances d < 3a the string fluctuations are not fully developped. The "Coulomb" 
behaviour of the total fleld energy observed in flg. ^ is mainly due to the parallel component of 
the magnetic field which gives the smallest contribution to the linear increase of the energy. This 
demonstrates that the behaviour of the flux tube cannot be explained easily within a classical 
static model. 



9 Conclusions 

Our motivation was a detailed quantitative analysis of the dual superconductor picture in U{1) 
lattice gauge theory. By considering the duality transformation of correlation functions we 
improved the understanding of the behaviour of expectation values of flelds and currents in the 
presence of a charge pair, as well as the significance of string fiuctuations. Moreover, simulating 
the dually transformed theory provided numerical flux tube data for low temperatures and large 
charge distances not obtained before in four-dimensional U (1) lattice gauge theory. 

We found that flnite temperature effects may be neglected for time extents exceeding the 
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Figure 8: The total field energy as a function of the charge distance d for NT=16 and NT=4. 
For higher temperature the periodic excitations in time direction mainly contribute to E'^ and 
therefore increase the total energy. 




Figure 9: The various contributions to the total energy (for NT=16). The parallel component 
of the electric field increases almost linearly for all distances, the parallel magnetic component 
develops its full slope after d = 3a only. 
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charge distance, while for smaller time extents the flux tube widens significantly. Further we ex- 
amined the electric field and magnetic current distribution and found that its description by the 
classical solution of Maxwell and London equation fails for large charge distances. The observed 
effects may be explained by fluctuations of the flux tube, on the other hand the current profile 
resembles that of a dual type-I superconductor. Roughening effects which become important for 
large distances are described well by the effective string picture of confinement which predicts a 
logarithmic increase of the fiux tube width. Finally, the analysis of the total field energy revealed 
the responsibility of the parallel magnetic field for the Coulomb behaviour at small distances. 

At this time, no final statement should be made which effective model is really able to 
describe strongly coupled U{1) gauge theory. Certain aspects of the flux tube may be explained 
by a "classical" dual superconductor picture, while others cannot. Further studies considering 
periodically closed flux tubes (torelons) are in progress and promise to answer some of the 
remaining questions. 
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